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Abstract. We investigate phase- insensitive linear amplification at the quantum limit 
for single- and two-mode states and show that there exists a broad class of non-Gaussian 
states whose nonclassicality survives even at an arbitrarily large gain. We identify the 
corresponding observable nonclassical effects and find that they include, remarkably, 
two- mode entanglement. The implications of our results for quantum cloning outside 
the Gaussian regime are also addressed. 
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1. Introduction 

When one amplifies a quantum state, a certain amount of noise is introduced as a 
fundamental requirement with profound implications for quantum communication and 
measurement; noise disallows the generation of perfect quantum clones pQ, which enables 
secure quantum communication, and it requires that quantum measurement — which 
necessarily involves amplification — be accompanied by decoherence, disallowing the 
mapping of a microscopic pure state to a macroscopic pure state. 

Following the work of Caves [2], which established quantum-noise limits for both 
phase-sensitive and -insensitive linear amplification, numerous studies have investigated 
how phase-insensitive linear amplification (PILA) washes out nonclassical effects. Hong 
et al. [3], in particular, showed that sub-Poissonian statistics and squeezing disappear 
below the quantum cloning limit (gain G — 2), and it has been argued that nonclassical 
effects disappear completely at sufficiently large gain [4]. The latter argument relies, 
however, on non-optimal amplification (e.g., by a bath of two-state atoms [5]); in 
quantum-limited amplification, states can remain nonclassical at arbitrarily large gain — 
e.g., this is known for Fock states [6] and Schrodinger-cat states [7]. It is known also that 
amplitude-squared squeezing [8] and certain higher-order photon statistics [9] can persist 
slightly above G — 2. More generally, which nonclassical states or nonclassical effects 
persist beyond the cloning limit has not yet been systematically studied. Moreover, the 
study of PILA has been confined mostly to single-mode states, and there are a number 
of crucial issues to address outside the single-mode regime; in particular, the question 
of whether entanglement, or nonclassical correlation, can survive is of broad interest to 
quantum information science [TO] . 

In addition to these consideration of PILA per se, linear amplification is generally 
the basis of quantum cloning schemes, particularly for continuous variables (CVs). A 
quantum doner based on PILA and a beam-splitter has been identified as the optimal 
Gaussian doner for coherent states and the scheme can be extended to other 
CV states as a general tool for quantum cloning. Little is known, however, about the 
behavior of non-Gaussian states under linear amplification. It is thus important to carry 
the study of PILA systematically beyond the treatment of the Gaussian regime found 
in earlier works [UJ [12] . 

In this paper we investigate both single- and two-mode states under the nonunitary 
(phase-insensitive) process of PILA; our treatment naturally extends to multi-mode 
states. As the unitary (phase-sensitive) amplifier is able to create nonclassicality, e.g., 
squeezing, it is ruled out in this work. We connect the critical gain at which all 
nonclassical effects disappear to the nonclassical depth defined by Lee [13] , and identify 
effects that remain observable at large gain. We show that there exists a broad class of 
single- and two-mode non-Gaussian states whose nonclassicality, remarkably including 
two- mode entanglement, survives even at arbitrarily large gain, in sharp contrast to the 
behavior known for Gaussian states [12]. This might also be contrasted with the recent 
work by Allegra et al. that presents evidence in support of the maximum robustness 
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of Gaussian entangled states against noisy thermal environments [3]. Furthermore, 
we show that the PILA-based quantum cloner can make more than two clones in a 
nonclassical state with fidelity better than classical scheme if the input has nonclassical 
depth larger than 1/2. Other aspects of PILA-based quantum cloning are also briefly 
addressed. 



2. Phase insensitive linear amplification 

We begin with some background to quantum-limited PILA. Consider a mode with 
creation and annihilation operators a) and a, which after amplification at gain G yields 

a G = y/Ga + VG^Tvl (1) 

In Eq. (pQ) v is a vacuum mode introduced to account for the fundamental noise in 
the amplification; do thereby becomes a legitimate field operator, with commutator 
[da, d G ] = 1. After taking the trace over the mode v , the action of the amplifier can be 
described by a completely positive map, with Kraus decomposition 

oo 

Pout = A[pin] = AnP^Al, (2) 
n=0 

where 

A n = (G - l) n/2 exp[-ln(G)ao72] (o + ) n / Vri. . (3) 

In practice, quantum-limited PILA can be achieved by injecting the signal mode a 
into a nondegenerate parametric amplifier (NDPA) whose idler is in the vacuum state. 
Remarkably, as in [151 US] > h can also be realized using only linear optics by employing 
homodyne detection followed by feed-forward. In this paper, the NDPA framework is 
mostly employed. Thus, the two-mode input to the amplifier is written p av = pi n £g> |0)(0|, 
where p in = / d 2 aP in (a)\a)(a\ is the signal- mode input to be amplified, with P; n (a) 
being its P-function. The amplified output is given by pc = tr v (UpAp av Up A ), where 
Upa = exp[Kt(a^v^ — av)} executes amplification with gain G = cosh 2 (/tt). After some 
calculation, the P-function of the amplified output is found to be 



P G (a) = ^ j d 2 f3P in (P) exp 



( (3-a/VG 2 



V 



1 - l/G 



(4) 



which amounts to a Gaussian convolution of the input plus rescaling of the phase space 
by y/G. 

Equation OH) establishes two important points about PILA: 
(i) In general, a single-mode state can be represented by an s-parametrized phase-space 
distribution, W s (a), with the parameter s in the range — 1 < s < 1 — e.g., Glauber- 
Sudarshan P-function (s = 1), Wigner function (s = 0), or Q-function (s = —1). This 
distribution is known to be well-behaved, with no singularity, for s < 0. In addition, 
two different distributions are related to each other by a convolution [17] : 



Linear amplification and quantum cloning for non- Gaussian continuous variables 4 



where s' < s. Therefore we see, comparing Eqs. (j3j) and (151), that the P-function always 
exists in a well-behaved form after amplification with G > 2. 

(ii) Equation (TJJ may be connected to the nonclassical depth defined by C. T. Lee [13J. 
The Glauber- P function becomes positive when smoothed by a Gaussian of sufficient 
width, 

P'(a) = ^pj d 2 (3P{f3) exp (~ \f3 - a] 2 ) . (6) 

The minimum r G [0, 1] required for a positive P-function is taken as a nonclassicality 
measure. Thus, by comparing Eqs. (jlj) and (|6]), we see that a state of nonclassical depth 
t remains nonclassical under PILA up to the critical gain G c — 1/(1 — r). This shows 
that Gaussian (squeezed) states, which have nonclassical depth in the range r < 1/2, 
cannot survive the cloning limit, i.e., for them G c < 2. 

In contrast to Gaussian states, "maximally" nonclassical states (r = 1) remain 
nonclassical at any gain, i.e., G c = 1/(1 — r) — > oo. It is therefore interesting to 
identify those states having r = 1. There exists a sufficient, though not necessary, 
criterion to detect maximally nonclassical states [18]: if p is orthogonal to a certain 
coherent state — i.e., if (a|p|a) = for some a — the state is maximally nonclassical. 
For instance, Liitkenhaus and Barnett [18] showed that all pure single-mode states, 
apart from Gaussian states, satisfy this condition. We thus conclude that all pure 
non-Gaussian states remain nonclassical under PILA at any gain G. Furthermore, any 
state (pure or mixed) can be transformed so that it becomes maximally nonclassical, 
and thus its nonclassicality persists to arbitrarily large gain; this is achieved by the 
photon- addition scheme, p' ~ a) pa [19]; p' necessarily satisfies (0|p'|0) = 0. 

2.1. Observable nonclassical effects 

In the above we have shown that there is a broad class of single-mode nonclassical states 
whose nonclassicality survives even at arbitrarily large finite gain. Let us now ask which 
nonclassical effects can be experimentally observed at high gain. If p is classical, its P- 
function is positive definite, so that we may construct, for any operator / = f(a\a), 
the positive quantity 

(:/*/:) = Jd 2 a\f(a*,a)\ 2 P(a)>0, (7) 

where : : denotes normal-ordering of the operators a^ and a. A violation of this inequality 
is an indicator of nonclassicality. We aim, then, to construct an operator T = pf that 
violates the inequality ([7]) for a state pa amplified at gain G. 

Since the P-function is well-behaved for G > 2, whenever pa is nonclassical at 
higher gain there must exist a finite region of phase-space (at least one), centered on a 
point «cb within which Pg((x) is negative. One can then readily construct a Gaussian 
function T(a) = exp(— \a — a \ 2 /cr), with variance a, such that J d 2 aPG(a)T(a) < 0. 
In operator form the negativity is expressed as 

(:f :) = {[(a - l)/ a ] (at - Q S )(a - ao) ) < 0, (8) 
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where we use the normal-ordering relation exp(8&b) = : exp[(e e — l)¥b] : for a field 
operator b [T7]. It serves as an observable nonclassical effect that persists to large gain. 
For simplicity, we may set «o = (equivalently, displace the state by — ao), in which 
case the defined nonclassical effect is a property of the photon-number distribution, 
P(n), and expressed by ([(a - l)/af & ) = E^° =0 [(^ - l)/o-] n P{n) < 0. It is testable in 
a photon-counting experiment. 

As an example, the one-photon Fock state, after amplification under PILA at gain 
G, has P-function 

which is negative for \a\ 2 < 1 — 1/G. Its nonclassicality is observed as ([(a— l)/cr] ata ) < 0, 
for any a < 1 and any G. In particular, by choosing a = 1/2 this becomes a condition 
on the photon-number parity — ((— l) ata ) < 0. 

Alternatively, one can construct a closely related nonclassical effect, more accessible 
in practice, that is insensitive to detector efficiency. We adopt a test operator in phase- 
insensitive form, / = Y,j CjCi^a? . A classical state must then satisfy Eq. (j7j) for arbitrary 
complex coefficients Cj, which yields the requirement DetfM*™)] > 0, for n = 1,2,... 
[91 [20] , where M^> is an n x n matrix whose elements, 

MW = (:(ota)^- 2 :) I (10) 

can be measured in photon-coincidence detection; Det[M^] < 0, for any n, indicates 
nonclassicality — e.g., with n = 2, Det[M( 2 )] < indicates sub-Poissonian photon counts. 

Photodetection with efficiency r] < 1 can be modeled by ideal detection after 
the signal mode, a, has been mixed with a vacuum mode, v, at a beam-splitter of 
transmissivity ^/rj. The detected mode, a d , is expressed by = ^Jrja + yT — rjv, which 
gives the matrix elements M-™\ et = r/ J+jf_2 ( : (a^a) J+J_2 : ). The determinant of this matrix 
turns out to depend on detection efficiency through an overall scale factor only — it is 
proportional to r] n2 ~ n — so the proposed nonclassical condition is formally (though not 
in the magnitude of the violation relative to experimental noises) insensitive to rj . One 
can show that the condition Det[M^] < holds up to a critical gain G^\ which can 
take a value far beyond G = 2. We plot G^ for an input one-photon Fock state in 
Fig. 1; the critical gain increases monotonically with n, which shows that higher-order 
nonclassicality survives a larger gain. 

3. PILA-based quantum cloning 

We now turn our attention to the cloning of continuous variable (CV) states using 
PILA. A quantum doner based on PILA followed by beam-splitting is identified as 
the optimal Gaussian doner for coherent states [TH [21] and has been experimentally 
implemented [TE] . It may not be possible to construct a universal optimal doner for the 
entire set of CV states, including all non-Gaussian states. Nevertheless, it does seem 
meaningful to take the PILA-based doner as a universal machine for a broader set than 
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Figure 1. Critical gain G c at which the nonclassicality condition Det[M^™^ < 
ceases to be valid for the one-photon Fock state under PILA. 



the Gaussian states alone. Suppose one amplifies a given input at gain G and injects 
the amplified output into a series of G — 1 beam splitters, each with transmissivity 
ti = 1/y/G + 1 — i, i = 1, . . . ,G — 1. In this configuration all the G output clones 
are reduced to an identical local state having the same mean amplitude as the input, 
(odone) = (Sin)- In such symmetric cloning, one readily sees that the P-function of each 
clone is directly obtained from that of the amplified state — phase-space scaling by 1/ VG 
due to beam-splitting — as P c ione(a) ~ Pc(VGa). That is, using Eq. (B|, 

^cione(«) = 1} / d 2 (3P m (P) exp ^- T _l_|/9 - «| 2 ) , (11) 

which, in view of Eq. (JSJ), is just the s-parametrized distribution of the input with 
s = 2/G — 1 > — 1. We now deduce the following: 

(i) In a classical measure-then-prepare strategy, where position and momentum are 
simultaneously measured to estimate an unknown state, the P-function at the output 
equals the Q-function at the input [22J. One can thus see that the PILA-based quantum 
doner is always superior, in terms of fidelity, to the classical scheme, regardless of 
the set of input states and number, G, of clones. The fidelity, P, between a pair 
of states, p\ and p 2 , is given by P = Tr{p 1 p 2 } = n' 1 f d 2 \Ci(\)C 2 (—\), where 
Cj(A) = Tr{Pj(A)pj} is the state characteristic function, with D^a) = exp(craj — a* hi), 
i — 1, 2, the displacement operator [T7]. When C 2 (X) is the same as the s-parametrized 
characteristic function of p\ — i.e., C 2 (A) = C*i(A) exp(s|A| 2 /2) — the fidelity is reduced 
to P = 7T _1 / <i 2 A|Ci(A)| 2 exp(s|A| 2 /2), which monotonically increases with s. Therefore, 
the fidelity between each clone and the input state decreases with increasing number, 
G, of clones, owing to the fact that s — 2/G — 1; it is, however, higher than that of the 
classical scheme with s — — 1. 

(ii) If the set of unknown input states has nonclassical depth r, all clones remain 
nonclassical up to the critical clone number G c = 1/(1 — r) [compare Eqs. (jSJ) and ( TTTT) ] . 
One cannot, for instance, have two clones of a squeezed state remain squeezed, since 
t < 1/2 does not give G c > 2. On the other hand, arbitrarily many clones of a state 
from the set of all pure non-Gaussian states can be created as nonclassical states since 
T = 1 =>• G c = oo. 
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4. Two-mode case 

We now extend our treatment to the two-mode case. In particular, we want to see 
how a nonclassical correlation, e.g., two-mode entanglement, evolves as each mode is 
amplified. An extension of Eq. (jJJ) gives the output two-mode P-function 

r 2 
P G (a u a 2 ) = I d 2 (3 l( l 2 (3 2 P in ((3i, fa) J[ Zi(a h A), (12) 

i=l 

where Zi = [n(Gi — l)] _1 exp —(1 — l/G.j) _1 |/3j — ai/\/G~i\ 2 , i — 1,2. One readily sees, 
from Eqs. ([5]) and ( |T2l) . that all Gaussian two- mode states, e.g., squeezed states, become 
classical for G > 2 as their original s-parametrized distributions are positive definite 
for s < [12] . Thus, if entanglement is to survive above the cloning gain, G = 2, its 
survival must be sought outside the Gaussian regime. 

For clarity we consider a particular but broad class of entangled states; namely, 
those obtained by combining a single-mode nonclassical state, p in , and a vacuum state 
at a beam splitter of transmissivity t = cos 9 [23] . This class serves as a useful tool to 
study many interesting entangled states: if |\l/i n ) is the Fock state \N), for example, one 



obtains a bimodal two- mode entangled state [24J, the simplest being the single-photon 
state cos#|l,0) + sin#|0, 1); if \^m) is the Schrodinger-cat state \a) — \ — a), one obtains 
the entangled coherent states |a;cos#, asin^) — | — a cos 9, — asin^). 

Let the input to the beam splitter be represented by the P-function P(ai,a2) = 
Pm{ a i)fi^ ( a 2) ■ Then, from the beam-splitter transformation, a[ = dicos0 + 5-2 sin 9, 
a! 2 = —Si sin^ + a<i cos 9, and Eq. ffl2l) . the amplified two- mode state has P-function 



where 



|/3-ai 2 /£i 2 | 2 ), (13) 

^ _ exp [\a 12 \ 2 /g l2 - K| 2 /(d - 1) - \a 2 \ 2 /(G 2 - 1)] 

12 " tt 2 (G 1 -1)(G 2 -1) 
Q 12 = \^[ 2 cos 2 9 + X^ 2 sin 2 9, 

1 1 

«i 2 = 7= oti cos 9 + a 2 sin 9, 



A^\/ G\ X 2 \G 



2 



with Aj = yl— l/Gi, i = 1,2. On the other hand, the local state of mode 1 has the 
P-function 



PgM = 1 J rf 2 /3P„(/3)exp| 



cos 2 9 



— ai/yGi cos 9 



2^ 



.(14) 



7r(Gi - 1) y V Af 

with cos^ — > sin9 in the similar expression for mode 2. Equations ( fl3|) and (fl4|) show 
how single- and two-mode nonclassicality decohere under PILA. Let us consider, in 
particular, symmetric PILA, i.e., G\ = G 2 = G. Then, if the input to the beam splitter, 
P-function P in (o;i), has nonclassical depth r, one sees from Eqs. ([6]) and f JT3|) that the 
global two- mode state remains nonclassical below the critical gain G c = 1/(1 — r). 
On the other hand, from Eqs. ([6]) and (fl4|) . the local states remain nonclassical up 
to G c i = 1/(1 — rcos 2 #) and G c2 = 1/(1 — rsin 2 ^), respectively. The local critical 
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values are smaller than the global critical value, G c i,2 < G c ; thus, for the range of gain 
satisfying G c i j < G < G c , the nonclassicality of the amplified state exists in the form 
of correlation exclusively. 

To take an example, when the entangled state \a, a) — | — a, —a) or |1, 0) + |0, 1) is 
amplified (r = 1), this analysis setting 9 = tt/4 shows that the nonclassical correlation 
survives at arbitrarily large gain, with G c = oo, whereas all local nonclassicality 
disappears above the cloning limit, G = 2. Since the amplified fields possess a well- 
behaved P-function at high gain, the method of [25] might be extended to the two-mode 
case to directly measure the P-function, thus demonstrating persistent nonclassical 
correlation for experimentally accessible non-Gaussian states. 

4-1. Persistent entanglement 

We must ask whether the surviving nonclassical correlation is indeed attributable 
to entanglement, since Marek et al. [26J find that separable states also can exhibit 
nonclassical correlation. It can be shown that the non-Gaussian states considered here 
are negative under partial transposition, a direct proof that they are entangled [27] . Let 
us illustrate persistent entanglement for the class of states generated by combining, 
at a 50:50 beam splitter, a maximally nonclassical single-mode state satisfying the 
condition (ct\p\ n \a) = [H] and a vacuum state. We prove that all such two-mode 
states remain entangled for any gain G. We may, without loss of generality, set a = 0; 
local displacement before the beam splitter is equivalent to local displacement after 
amplification, which does not change the entanglement. Thus, the input-state condition 
can be taken as (0|p; n |0) = 0, i.e., vanishing vacuum-state probability. We denote by iV 
the lowest nonzero occupation number of p in . 

Note now that the actions of the beam-splitter and PILA commute in symmetric 
amplification, as one readily verifies in the Heisenberg picture. Furthermore, the 
condition (j\p- m \j) = 0, j < N, holds also for the amplified state, Pi n ,G) which can 
be confirmed from Eq. (T5]). On the other hand, after amplification at gain G a vacuum 
state becomes a thermal state, |0)(0|c = G -1 X)n(l — 1/G) n \n)(n\. Therefore, for the 
two-mode amplified state, pe,g = Ubs • pm,G ® |0)(0|g • U^ s , we obtain 



Equation ( !T5|) gives a zero diagonal element (0, 0|/9e,g|0, 0) = and nonzero off- 
diagonal element (0, N\p Ej c\N, 0) ^ 0. The entanglement of pe,g is then proved by 
negativity under partial transposition in the subspace spanned by |0, 0) and \N, N); the 
determinant of pe,g under partial transposition is (0, 0\pe,g\®, 0)(iV, N\pe,g\N, N) — 



{0J\pE, G \j,0) = ^(0,0 
7 




(15) 



(0,iV|p£, G |iV,0> 2 <0. 
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To further verify this persistent entanglement, known methods based on partial 
transposition can be used [28] — e.g., for amplification of the state |1, 0) + 10, 1), a witness 
operator W = |e_)(e_| PT , with |e_) = |0,0) - (v / TTA T - A 2 )|l, 1) and A = y/l - l/G, 
is constructed to show that Ty{Wpe,g} < 0. 

4-2. Cloning of quantum correlation 

Finally, we briefly discuss the case where PILA followed by beam splitting is applied as 
in Sec. Ill to each of two modes, to achieve quantum cloning of nonclassical correlation 
or entanglement. Our analysis shows that if the PILA-based quantum doner makes 
multiple copies of a two-mode entangled state, any pair of clones possesses a nonclassical 
correlation up to the critical clone number G c — 1/(1 — r); the beam splitter following 
PILA simply rescales the P-function, leaving its negativity unchanged. On the other 
hand, negativity under partial transposition disappears at the two-clone limit (G = 2) 
for symmetric cloning, due to vacuum noise at the beam-splitter. A subtle contrast 
between nonclassical correlation and entanglement therefore emerges in the symmetric 
quantum cloning of two-mode non-Gaussian entangled states [26] . 

If, however, only one party makes clones (G2 = 0), it turns out that entanglement 
survives between any clone of mode 1 and the original copy of mode 2 regardless of G\ 
for both Gaussian and non-Gaussian pure inputs [29]. This can be explained in a general 
framework of one-sided quantum channels for CVs and more details will be investigated 
elsewhere. 

5. Concluding remarks 

We have shown that certain nonclassical effects persist far beyond the cloning limit 
(at gain G = 2) in quantum-limited PILA. We identified a broad class of single-mode 
and two-mode non-Gaussian states that survive arbitrarily large gain and associated 
observable effects that are testable by photon counting. Remarkably, these persistent 
nonclassical effects include quantum entanglement for two-mode non-Gaussian fields. 
This result is to be contrasted with that of Ref. [H], which provided some evidence 
for the maximal robustness of Gaussian entanglement at fixed energy in the presence 
of thermal noise. Although it was further conjectured that Gaussian entangled states 
would generally be maximally robust under noisy Markovian interactions, our work 
provides a counter- argument; Gaussian entanglement disappears above G=2, whereas 
non-Gaussian entanglement persists beyond G = 2 under the Markovian evolution of 
PILA [5]. Our results also show that a quantum doner based on PILA can produce 
many nonclassical clones with better fidelity than the classical measure-then-prepare 
scheme for a set of non-Gaussian states. Demonstration of these results seems feasible 
in view of recent experimental progress [El [161 EE] . 
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